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ABSTRACT 

We present deep MUSE integral-field unit (IFU) spectroscopic observations of the ^ 150kpc 
Lya halo around the z = 4.1 radio galaxy TNJ J1338 —1942. This 9-hr observation maps the 
two-dimensional kinematics of the Lya emission across the halo. We identify two Hi ab¬ 
sorbers which are seen against the Lya emission, both of which cover the full ^ 150 x 80kpc 
extent of the halo and so have covering fractions ~ 1. The stronger and more blue-shifted 
absorber (Av ~ —1200kms“^ from the systemic) has dynamics that mirror that of the un¬ 
derlying halo emission and we suggest that this high column material (n(Hl) ~ 10^®'^ cm“^), 
which is also seen in CiV absorption, represents an out-flowing shell that has been driven by 
the AGN or the star formation within the galaxy. The weaker (n(Hl) ~ 10^^ cm“^) and less 
blue shifted (Av ~ —500kms“^) absorber most likely represents material in the cavity be¬ 
tween the out-flowing shell and the Lya halo. We estimate that the mass in the shell must be 
of order 10^° M© - a significant fraction of the ISM from a galaxy at 2 ; = 4. The large scales 
of these coherent structures illustrate the potentially powerful influence of AGN feedback 
on the distribution and energetics of material in their surroundings. Indeed, the discovery of 
high-velocity (~ 1000km s“^), group-halo-scale (i.e. > 150kpc) and mass-loaded winds in 
the vicinity of the central radio source are broadly in agreement with the requirements of 
models that invoke AGN-driven outflows to regulate star formation and black-hole growth in 
massive galaxies. 

Key words: galaxies: evolution — galaxies: high-redshift — radio 


1 INTRODUCTION 

Theoretical models of galaxy formation are increasingly suc¬ 
cessful in their ability to predict the observed properties of 
galaxies across cosmic time. The most successful models in¬ 
voke a dramatic form of energy injection from the central ac¬ 
tive galactic nuclei (AGN) and from star-formation and super¬ 
novae to drive high-velocity (lOO’s kms“^), galaxy-wide (i.e. 

lOkpc) and mass-loaded energetic outflows. These winds ex¬ 
pel significant amounts of gas from the host galaxy and con¬ 
sequently this shuts down future star formation and black hole 
growth. This process should also enrich the larger scale en- 


vironment with metals (e.g. Silk & ReesI 19981: FabianI 

1999; 

Benson et alJ 2003: 

Granato et al. 2004: Di Matteo et all 

2005; 

Snrinsel & Hernauist 

12003: IScrinsel et alJ 2005: 

Honkins et al. 

20061 20081: Booth & Schavel 2010l: iDebuhr et alj 

2012h. 

In the 


models, this feedback mechanism is required to prevent the run¬ 
away cooling of gas in mass ive halos and so match the galaxy lumi - 
nosity function at 2 = 0 (e.g. [Bower et al.l2006LICroton et al.l200^ . 


If this feedback mechanism is the dominant route by which 
massive galaxies terminate their star formation, then it is critical 
that the AGN feedback is at its maximum efficiency at the epoch 
where the bulk of the stars in th e most luminous el liptical galaxies 
appear to be formed, 2 > 3-5 dNelan et al.ll2005h . This is there- 
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fore the era to search for evidence of AGN feedback in massive 
galaxies. High-redshift radio galaxies are one potential laboratory 
to search for evidence of this feedback “in action”, as they appear to 
be amongst the most evolved and most massive galaxies in the early 
Universe, host energetic AGN, and have the prope rties expected for 
the pr ogenitors of today’s central cluster galaxies dPe Breuck et al.l 

I2OOA . 

Potential evidence for such “AGN feedback” may come from 
the discovery of very extended Lya halos with spatial extents 
of ~ 100-200 kpc which are being found in increasing numbers 
in the high-r edshift Universe, including around high-red shift ra¬ 
dio galaxies dvan Oiik et al.l 1 19971 : IVenemans et alJl20()% . Many 
of these appear to be associated with young, massive galaxies 
which display high-levels of starburst an d AGN activity (e.g. 
Ivan Breugel et al.ll998l : lMatsuda et al.l201lh . One simple interpre¬ 
tation is that these Lya halos represent gas reservoirs which are 
being ionised by UV emission, shocks and outflows from the active 
galax i es within the halos (e.g. iTaniguchi et al.|[200ll : lGeach et al.l 
I 2 OO 9 I : ISmail & Blundelll [20131) . However, theoretical models have 
also suggested that these halos may instead b e powered by cool- 
from gas in so-called “cold-flows” (e.g. Ipardal et al.ll200ll : 




IPekel et aUbOOgh . The spatially resolved dynamics and line pro¬ 
files of these Lya halos, as well as their chemical enrichment 
should be able t o differentiate which mechanism is respons ible for 
their formation dSower et aDl2004 ISmail & Blundeiill2013h . 

The t wo-dimensional Lya surface brightness distribu¬ 
tions (e.g. iMatsuda et al.l 12012h or one- dimensional kinemat- 


Binette et alj2000l;IVillar-Martm et alj2003l;l0h' 

/ama & 

Taniguchi 

2004 

Matsuda et al.l2006l; Humnhrev et al.l2007 

. 20091; 

Yang et al. 


in shaping the environments around massive, high redshift galax¬ 
ies on <: lO’s of kilo-parsec scales. For example, observations have 
reveals that the spatial extent of the Lya-emitting gas is correlated 
with the spatial scales of the jets, whilst the distortion seen in the 
Lya morphology is also typically correlated with the distortion in 
the radio emission (the extended emission line region morphol¬ 
ogy always appears to be b rightest on the side of the radio lobes 
that i s closer to the nucleus: [McCarthy et alJl9^ : lTadhunter et al.l 
l200(]h . suggesting that powerful radio sources are (at least partly) 
responsible for environment al effects on large sc ales. In terms of 
the large scale environment, Ivan Oiik et alj ( Il997t) showed that the 
galaxies with radio jets that are extended on scales smaller than 
^ 50 kpc have the strongest strong Hi absorption seen against Lya 
(and likewise, those with radio extents ^ 50 kpc do not show strong 
absorbers). One explanation is that when radio sources growth they 
destroy the Hi. However, this model also requires a duty cycle 
in which the “feedback” expels gas, which cools then to form a 
shell of Hi. a subsequent “feedback” event then ionises the gas. 
This model may also explain why the dynamics of the gas in 
the inner halo (located within the radio source) tend to be turbu¬ 
lent (cr ^ 700-15 00 km s~^) whilst th e out halo is more quiescent 
(a ~300kms~^: lvan Oiik et aDll996l) . 

To measure the energetics of the Lya emission, search for 
and place dynamical constraints on the out-flowing material re¬ 
quires resolved spectroscopy. However, relatively few galaxies 
have been studied using their two-dimensional Lya velocity field, 
as required to truly understand their origin. The first attempt to 
derive the mor phologies and dynamics of exten ded Lya halos 
were made by iHippelein & Meisenheimeil ( 119931) who found a 
coherent and spatially extended (^lOOkpc^) and moderate col¬ 
umn, n(Hl)~10^® cm^, Hi absorber in front of a 4C radio galaxy 


at 2 = 3.8 (see also lAdam et al.l [l997l) . However, relatively few 
galaxies have been studied using their two-dimensional Lya ve- 
locity field, as required to truly understand their origin (e.g. see 
Boweret^ 2004; Wdman et alj 2005j: Villar-Martin et al. ^074 


Bgnchez & Humphrey 20091 : Weiimans et al. l20ld: iMartin et*^ 

2014 for a few examples). In part, this is due to the high red- 


shift and hence large luminosity distances required to redshift Lya 
into the optical band, which means that spatially resolving the 
low surface brightness features in the outer halo have, until now, 
been beyond the sensitivity limits of most instruments (even on 8- 
meter telescopes). However, with the commissioning of the sensi¬ 
tive Multi-Unit Spectroscopic Explorer (MUSE) on the ESO / VLT, 
it is now possible to carry out these observations and so determine 
the nature of extended Lya halos. 

In this paper, we exploit MUSE observations to determine 
the nature of extended Lya halo around a powerful radio galaxy, 
TNJJ1338-1942 at 2 = 4.1 (hereafter T NJJ1338). TNJ1338 wa s 
identified as hosting a giant Lya halo by Venemans etadjjgOOd). 
In a s ubsequent paper. [^nemans et alj(l2007l) (see also Milev et d] 
|2004 identified twenty Lya emitters within a projected dis¬ 
tance of 1.3 Mpc and 600kms“^ of the radio galaxy (an over¬ 
density of 4—15 compared to the average field at this redshift), 
suggesting a mass for the structure of M~5 x lO^'^M© which 
makes structure a likely ancestor of a massive galaxy cluster at 
the present day (with the radio galaxy likely to be a brightest 
cluster galaxy in formation). This high-redshift radio galaxy re¬ 
mains one of the the brightest and most luminous known in Lya 
{Lhya ~ 5 X 10^^ erg s~^). Both its Lya profile and radio structur e 
are very asymmetric jPe Breuck et al.iri999l : Iwilman et aDl2004h . 
which indicates strong interaction with dense gas, and the rest- 
frame radio luminosity is comparable to that of the most lumi¬ 
nous 3C radio sources (Pizsmhz ~4x 10®® erg s“^ Hz“^ sr~^). 
The MUSE observations presented here allow us to map the emis¬ 
sion within the ~ 150-kpc long Lya halo, matching the sensitivity 
of previous narrow-band Lya imaging of this structure, but with 
the crucial advantage of high spectral resolution kine matics, neces¬ 
sary to differentiate the signatures of infall or outflow jAdams et al.l 
l2009l : lBower et alj|2004rwilman et alj2005l) . 

The structure of this paper is as follows: in § 2 we present 
the MUSE observations and data-reduction. In § 3 we discuss the 
morphology and energetics of the Lya halo and those of the Hi ab¬ 
sorbers seen against the underlying emission. In §4 we present our 
conclusions. Throughout this paper, we adopt a cosmology with 
flm =0.27, Da =0.73 and Hq =71 kms“^ Mpc“® giving an angu¬ 
lar scale of 7.0 kpc arcsec~^ at 2 = 4.L 


2 OBSERVATIONS AND DATA REDUCTION 

We observed TNJJ1338 with MUSE for a total of 9.0hours be¬ 
tween 2014 April 30 and 2014 June 30 during a combination of 
commission and Science Verification time. The MUSE IFU pro¬ 
vides full spectral coverage spanning 4770-9300A, a contiguous 
field of view of 60 " x 60 ", with a spatial sampling of 0.2” pixel“^ 
and a spectral resolution of R = 3000 at A = 6200A (the wavelength 
of the Lya at 2=4.1). The observations were carried out in dark 
time, good sky transparency and 0.7” U-band seeing. Each 1 hr ob¬ 
servation was split in to 1800 second exposures with a small (2”) 
dither. We obtained two overlapping pointings; the first centered at 
a: 13 38 28.0 5: —1942 24 and a second, offset by approximately 
30” at a: 13 38 25.0 <5: —1942 32. Both pointings cover the 2 = 4.1 
radio galaxy and the final field of view is L5'x 1.0' (Fig.[TJ. 


(c) 0000 RAS, MNRAS 000, 000-000 


































































































MUSE observations of TN J1338—1942 at z = 4.1 3 


cn 

"D 

c 

o 

o 

CD 

CO 

o 

< 


30 

20 

10 

0 


-10 


-20 


-30 


-40 



45 30 15 0 -15 -30 

Arc Seconds 

Figure 1. Colour, (V, I, LyQ 2 = 4 .i) image of the field around TNJ J1338 from the MUSE observations. This colour image illustrates the full extent the giant 
halo around the 2 : =4.1 radio galaxy. Moreover, a search for other emission lines within the MUSE cube identifies a new giant (~60kpc) halo around a 
foreground z = 3.31 galaxy located to the south west of the radio source. We also mark the positions of the two z = 4.1 galaxies in the halo of the radio galaxy 
(open diamonds). The origin of the image is centered on the radio galaxy at a: 13:38:26.1, <5: —19:42:30.5 (J2000). 


To reduce the data, we use the MUSE pipeline through ES- 
OREX which wavelength calibrates, fl at-fields the spectra and forms 
the datacuhe dWeilhacher et alj2014l) . Due to the excellent fiat-field 
stability of MUSE, sky subtraction was performed on each indi¬ 
vidual cube by identifying and subtracting the sky emission using 
blank areas of sky at each wavelength slice (i.e. no additional sky- 
only observations are required), and the final mosaic was then con¬ 
structed using an average with a 3-a clip to reject cosmic rays af¬ 
ter registering the stars in continuum images from each individual 
cube. 

A visual inspection of the final datacuhe reveals a very strong 
detection of the extended Lya emission from the z = 4.1 halo at 
6200A as well as Lya emission from a number of galaxies at the 
same redshift. In Eig. [T]we show a colour image of the field gen¬ 
erated from the MUSE cube. In this image, we use the continuum 
from the wavelength ranges 4900-6900A and 6900-9000A as the 
blue and green images respectively, whilst the red image is gener¬ 
ated by collapsing the cube over the wavelength range 6I00-6200A 
(which covers the Lya from the z = 4. 1 radio galaxy). The strong, 
extended Lya emission can clearly be seen, with a spatial extent of 
~ 150 X SOkpc. 

We note that a search for line emitters within the cube also 


detects a second strong, spatially-extended Lya halo around a red 
2 = 3.31 galaxy located approximately ~30" to the south west of 
TNJJI338 (Fig.E) which appears to lie in an over-density com¬ 
prising a number of other star-forming galaxies and AGN at this 
redshift. In addition, a large number of foreground [Oil], [Olll] and 
Lya emitters are detected in this volume. These will be discussed 
in a subsequent paper (Urrutia et al. in prep). 


3 ANALYSIS & DISCUSSION 

3.1 Morphology and Energetics of the Lya Halo 

The morphology of the Lya emission (Fig. [Til is very similar to 
that se en in the FORS 2 narrow-band imaging from IVenemans et ahl 
(l2007h . which has similar surface brightness limits as the Lya 
narrow-band image created from the MUSE cube (collapsing the 
cube over 120A around the Lya yields an image with a 5-a line 
flux detection in a 1.2" aperture of /so- = x 10~^® erg s”*^ cm“^ 
which is similar to that of the narrow-band image from 
IVenemans et al.ll2007l : = 6 x 10“^® ergs”*^ cm“^). Indeed, the 

morphology of the source shows a strong, spatially extended core 
which has a diameter of approximately 5” with a lower surface 
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Figure 3. Upper panel: Continuum subtracted, wavelength collapsed, narrow-band image of the TNJJ1338 Lyo halo from the MUSE cube (collapsed over 
the wavelength range 6100-6300A). The cross de note the position of the U V-continuum (and radio) counterpaifs, and the open triangles show the positions 
of the radio hot-spots (end points of the jets) from lPe Breuck et alj jl999ll indicating the compact size of the radio source compared to the highly-extended 
Lya halo. We also show example spectra of regions within the halo. These are typically extracted from 2” X 2” regions, and in all cases, the line profile is 
best-fit by two H-I absorbers superimposed on an underlying Gaussian emission profile. Both absorbers have a covering factor of unity; the first absorber, 
Hl-B appears blue-shifted from the (underlying) Gaussian emission line profile by ~ —700 km s~^ and has a column density of n(Hl)~ 10^®'® cm~^ whilst 
the second absorber, Hl-R is much closer to the underlying halo emission. An ~ —200 km s~^ but has a much lower column density, n(Hl)~ 10^^ ® cm~^. 


brightness halo which is extended along a roughly North-South di¬ 
rection which extends over a ~ 20” (150 kpc) in extent. The elon¬ 
gated Lya morphology is also consistent with the direction of the 
radio jets (Fig.[TJ, although the strongest jet emission is contained 
within the central 10 kpc, which may be indicative of previous 
AGN activity (which ionised the extended halo and which is now 
seen in Lya emission as the gas cools). We note that the elongated 
morphology of the Lya - in the same direction as the radio emis¬ 
sion - is consis tent with that seen in a number of other high-redshift 
radio galaxies tebttgering & Pentericci|[l99^ . 

As an order-of-magnitude estimate of the total Hll mass re¬ 
sponsible for the Lya emitting halo, we can assume case B recom- 
bin ation and that the gas is at a temperature of 10^ K, and follow¬ 
ing iDe^BreuckeT^ ( Il999h . using the total Lya emission intensity 
(corrected for absorption, see S I3.3b . the total Hll mass is approx¬ 


imately M(Hll)~ 10^ (/5 L 44 V 7 o)°'® M( 7 ), w here /s is the volume 
filling factor of 10“® llMcCarthv et alJfl9^ . Z /44 is the Lya lumi¬ 
nosity in units of 10 ^“* erg s“^ and V 70 is the volume of the region in 
units of 10^^° cm®. Adopting a size of 150 x 80 x 80kpc we derive 
M(Hll)-4 X 10® Mq. 

The spatial extent of the Lya halo in TNJJ1338 is compa- 
rabl e to that seen in a number of other high-re dshift radio galax¬ 
ies Jvan Oiik et al.lll99^ ; lHumDhrev et alj|2007l) . Indeed, extended 
Lya halos are common around high-redshift radio galaxies, ex¬ 
tending up to ~ 250 kpc - fa r beyond the continuum s tructures re¬ 
vealed by HST images (e.g. Ivillar-Martin et al.|l2003h . Moreover, 
the spatial extent of the Lya emitting halo around high-redshift 
radio gala xies also appear to c orrelate with the size of the radio 
emission jvan Oiik et al. I ll997h (although of course, this is likely 
to be a selection effect since the larger halos are likely to be more 
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Figure 2. The spectrum extracted from the core of the Lya halo compared 
to that of the extended halo. The core spectrum shows strong emission lines 
associated with the AGN activity (e.g. Nv, [OlV], CiV and Hell). The emis¬ 
sion from the halo (which should exclude any seeing-related contamination 
from the core, but will not exclude scattered emission within the halo) also 
shows weak CiV emission. The Lyo image of the halo (inset) shows the two 
regions from which the core and halo were extracted (for the core, all pixels 
within the inner contour were included and for the halo, all pixels between 
the second and third contours were included). The second inset compares 
the line profile of the Lyo, CiV and Hell from the core region. In this panel, 
we use the Hell emission to define the systemic redshift, 21 = 4.1057 (and 
transform the Lyo, CiV and Hell to the same rest-frame velocity scale). In 
the core region, the Lyo and CiV share the same characteristic (asymmet¬ 
ric) line profile, implying the presence of CiV in the absorbing material and 
suggesting that any absorbing materials must be enriched in metals 


easily detectable around powerful active galaxies since the active 
nucleus provides a more intense supply of ionising photons). Blind 
narrow-hand surveys have also identified giant Lyo emitting ha¬ 
los aro und star-forming gala xies with (relatively) low luminosity 
AGNs jMatsuda et aljEoTlI) as well as in structures c ontaining 
multiple radio-quiet quasars (e.g. ICantaluno et aflEoidh . suggest¬ 
ing that while an AGN may play an important role in providing the 
ionising source of radiation to ionise the halo, it is not a require¬ 
ment. Indeed, the Lya blohs (LABs) also have spatial extents up 
to ISOkpc and the Lya emission can he attributed to either pho- 
toionisation from t he AGN or star-formation from the host galaxies 
jGeach et alj200(^ . 

What powers the Lya emission in the extended halo of 
TNJ 1338 is unclear. As Fig. [T] shows, the Lya halo extends over 
lOx that of the radio jets detectable at GHz radio frequen¬ 
cies, which are contained within the central lOkpc. One pos¬ 
sible interpretation is that the emission in the extended halo is 

S lowered by the cooling of pristine gas within a dark matter halo 
Fardal et aDl200ll : iNilsson et alJbOOhl) . However, the far-infrared 
luminosity of TNJ 1338 implied from Herschel / SPIRE and 
JCMT/ SCUBA-850^m observations is Lir. =7 ± 1 x 10^^ Lq 
(SFR = 700Mq yr“^; lOrouart et al.l l2014h which gives a ratio 
of total Lya/IR luminosity of ~1 x 10“*® /4 x 10“*® =0.03. This 
luminosity ratio is similar to that from the A GN luminos¬ 
ity (Lagn ~2x 10^®ergs~^: IProuart et al.ll2014l) which gives 
Lnya/L agn = 0.05. Thus, only a small fraction of either the in¬ 
frared (star-formation) or AGN luminosity is required to power 
the Lya emission. However, we note that the high Lya/Hen ra¬ 


tio averaged over the halo is Lya/HeIlA1640= 18.7 ±0.5 (see 
SI3.2b which is difficult to e xplain with AGN photo-ionisation only 
jVillar-Martm et ^l2007al) . and suggests we are likely to be ob¬ 
serving a combination of energy sources that are heating the halo. 


3.2 Spatially Resolved Structure 

Before determining the velocity structure of the halo, we first de¬ 
fine the systemic redshift of the galaxy by extracting a spectrum 
from a 1 X 1” region centered on the UV (and radio) continuum 
counterpart. Unfortunately, no strong photospheric absorption fea¬ 
tures are present, but there are strong detections of Lya, CivA 1549 
and Hell A 1640 (see Fig. O. Since the Lya and CiV are reso¬ 
nant lines, we use the Hell emission line and fit a single Gaussian 
profile, deriving a redshift of 2 = 4.1057 ±0.0004 and a FWHM 
of 35 ± 3 A (or 1250 ± 100km s“^). In all following sections, we 
adopt 2 = 4.1057 as the systemic redshift when comparing veloci¬ 
ties or velocity differences. 

To investigate the spectral properties of the core and halo, 
we begin by extracting two spectra, one from a circular 3 x 3” 
region centered on the continuum counterpart (the core of the 
galaxy), and a second from the extended halo (excluding the core 
region), and we show these two spectra in Fig.[2 The core spectrum 
shows strong emission lines associated with the AGN (e.g. Nv, 
[OlV], Civ and Hell). The emission from the halo (which should 
exclude any seeing-related contamination from the core, but will 
not exclude scattered emission within the halo) also shows weak 
Civ emission. The detection of heavy elements (in particular CiV) 
in the ex tended halo suggest s that the halo is already m etal-enriched 
(see also lBinette et alj2006l : IVillar-Martm et alj2003l) . 

The line ratios from Carbon, Nitrogen or Helium compared 
to Hydrogen can be used to infer the physical conditions within 
the diffuse halo. Although the CiV and Hell emission lines lie in 
regions of high atmospheric OH emission, we derive line ratios 
(or limits) of [CiV] / Lya = 0.082 ± 0.003 and 0.064 ±0.013 for the 
“core” and “halo” respectively. For the Hell/Lya ratio, we derive 
Hell/Lya = 0.044 ±0.003 and <0.075 for the “core” and “halo” 
respectively - consistent with stellar or power-law ph oto-ionisation 
models with metallicities ranging from 0.05-1 Z/Z© dBinette et al.l 
l2006l:IVillar-Martm et alj2007al) . 

Next we extract a number of spectra across the halo by in¬ 
tegrating 2" X 2" regions and show these in Fig. Whilst the 
strong Lya emission is clearly evident, it is also apparent that the 
emission line is asymmetric. One possible interpretation of the Lya 
emission line profile is that it naturally arises from the transfer of 
Lya radiation through a neutral screen. If the Lya photons have 
to travel through a region of neutral hydrogen, they will resonantly 
scatter until they Doppler shift into the wings where the opacity is 
lower. Such a phenomenon will give rise to either a double-peaked 
profile dAdam jl975l : [Diikstra etalj 20061) . or, depending on the ge¬ 
ometry and kinematics of the neutral gas, an asymmetric profile in 
which one of the peaks is stronger. However, for this to be a viable 
explanation of the observed line profile in TNJ 1338, we require 
the Lya to be embedded in a (partial) cocoon of neutral hydro¬ 
gen of large optical depth. On the other hand, the asymmetry of 
the emission line, and the sharpness of the deep trough(s) suggest 
that Hi absorption is far more probable (e .g. iRottgering et al.l 19951 : 
Ivan Oiik et alJl9^ : IWilman et al.ll200^ . 

Within the central regions of the halo, the asymmetric Lya 
pro file was first studied using high-resolution UVES spectroscopy 
by dWilman et alj|2003) who showed that the line emission is best 
fit with an underlying Gaussian profile with two Hi absorbers with 
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Figure 4. Top row: Intensity map of the underlying Lya Gaussian emission line profile (Left), and spatial column density distribution of Hl-B (Middle) and 
Hl-R (Right). Bottom Row: Velocity structure of the Gaussian profile, Hf-B and Hl-R. In all cases, the velocity scale is set by the systemic, as measured from 
the Hell emission: z =4.1057. The solid point in the top left corner of each panel denotes the 0.7" FWHM seeing disk. In each panel, we also show the position 
of the radio core (and continuum counterpart of the galaxy; cross) and the radio lobes (open triangle). The underlying Gaussian emission line profile shows a 
maximum velocity gradient of ~ 1000 km s “ ^ across the 150 kpc (in projection). At the position of the galaxy, the Lya is blue-shifted by ~ 150 km s “ ^ from 
the systemic. The velocity dispersion profile of the underlying emission (which has an average value of tr = 570 ± 13 km s“^) is peaked at the position of the 
radio galaxy and the radio-jets, which may reflect the interaction of the radio plasma and the gas. Both Hf-B and Hl-R have covering factors of unity. Hl-B 
appears blue-shifted from the (underlying) Gaussian emission line profile by ~ —700kms“^ and has an average column density of n(Hl)~ 10^®’® cm~^. 
This material is most likely associated with an out-flowing shell of material, and we estimate the mass in the shell must be of order ~10^*^ M©- The second 
absorber, Hf-R is much closer to the underlying halo emission in velocity space. An ~ —200 km s~^ but has a much lower column density, with an average 
of ?i(Hl)~ ® cm“^. This material is most likely to be associated with the low-density ICM between the out-flowing shell and the Lya halo. 


column densities n(Hl) ~ 10^® ° cm“^ and n(Hl) ~ 10^®'^ cm~^ 
for Hl-B and Hl-R respectively. As we shall see, a similar com¬ 
bination of absorbers is required to reproduce our spectra across all 
of the halo. 

To describe the line profile, we begin by fitting a Gaussian 
profile emission line with a Hi absorber (which is described by 
a Voigt-profile). However, in all of the regions shown in Fig. 
this fit fails to provide an adequate description of the line, and in¬ 
stead two Hi ab sorbers are required to adequately describe the line 
profile (see also lWilman et al.|l2004h . The first of these absorbers, 
which we label Hl-R (the most redshifted) is typically offset from 
the central (underlying) Gaussian profile by —270kms“^, has 

a column density of rifHO-^ 10^^ ® cm~^ and Doppler parame¬ 
ter of &=100±4kms~^ whilst the second absorber, Hl-B (the 
most blue-shifted) is typically offset from the underlying Gaus¬ 
sian profile by 1200kms“^ but with a much higher column 
density, n(Hl)~ 10^®’® cm“^ and 6= 170 ±5 km s“^. Within the 
halo, the underlying Gaussian profile has an average velocity dis¬ 


persion of (T = 570 ± 13 km s~^, in agreement with the velocity dis¬ 
persions of the Hell emission line, and much broader than the two 
absorbers. The velocity dispersion of the underlying emission is 
peaked at the position of the radio galaxy (and the radio-jets), with 
cr = 750 ± 20 km s~^. This may reflect two components of the halo: 
the inner halo where there is a strong interaction between the ra¬ 
dio plasma and the gas, and th e outer (more quiesc ent) halo, lo¬ 
cated outside the radio source (Ivan Oiik et al. I l 19961) . In contrast, 
the much smaller velocity dispersion of the absorbing gas is more 
readily explained if these absorbers lie at relatively large distances 
from the parent galaxy. We will discuss the spatially resolved prop¬ 
erties of these absorbers in 8) 13.31 

Finally, we note that Ivan Oiik et al.] d 19971) showed that most 
high redshift radio sources with radio jets within ~ 50 kpc of the 
central galaxies show strong Hi absorption against their Lyo? emis¬ 
sion. This may suggest that when the radio sources grow, they de¬ 
stroy the Hi. Since the radio emission in TNJ 1338 is located within 
10 kpc of the galaxy, this radio source clearly falls into first cat- 


(c) 0000 RAS, MNRAS 000, 000-000 
























MUSE observations of TN J1338—1942 at z = 4.1 1 


egory. Within this scenario, there must be a duty cycle with the 
feedback expelling Lya-emitting gas, which cools to form a shell 
of Hi around it, and a subsequent feedback event ionises the gas. 


3.3 Dynamics of the Hi Absorbers 

As can be seen from Fig. one of the striking results from 
our MUSE observations is the coherence and spatial extent of 
both absorbers. The most blue-shifted absorber, Hl-B, which has 
n(Hl)~ 10^®'® cm~^ appears to be spatially extended across the 
entire 150 x SOkpc structure. The second absorber, Hl-R, which 
has a significantly lower column density, n(Hl)~ 10^'* ® cm“^ also 
appears to be extended over the same area. Thus, the cover¬ 
ing factor of both absorbers appears to be close to unity. Evi¬ 
dence for extended reservoirs of neutral gas have been seen in 
a number of giant Lya halos around high-redshift radio galaxies 


I van^)jik_e^lj [lj97l : iBinette et al.ll2000l : fvillar-Martm et al J2003l: 

Wilnian_et^|200 ^ as well as star-forming gal axies iBower et al.l 

20041 : Wilman et al.ll200^ : IWeiimans et al.ll2010l) . 

To examine whether the Hi absorbers are “pristine” - as sug¬ 
gested by the cold-flow model, in Eig.|2]we compare the Lya, CiV 
and Hell emission line profiles (transformed to the same rest-frame 
velocity scale). Although the CiV emission is lower signal-to-noise, 
it is clear that the Lya and CiV share the same characteristic (asym¬ 
metric) line profile, implying the presence of CiV in the absorbing 
material (without other species, it is difficult to derive the ionisa¬ 
tion of the absorbing material, but it appears that this absorbing gas 
must be metal enriched, arguing against a model in which the Lya 
results from the cooling of pristine material). 

With the two dimensional data provided by MUSE, next, we 
construct a spatially resolved map of the dynamics and column den¬ 
sity across the halo. To achieve this, at each spatial pixel, we av¬ 
erage the surrounding 1" x 1" region, and measure the signal-to- 
noise (S /N). If the S /N> 50, we then we attempt to fit the emis¬ 
sion line profile with an underlying Gaussian profile with two Hi 
absorbers, allowing the centroid, intensity and line-width (FWHM 
in the case of the Gaussian profile and Doppler 6-parameter in the 
Voigt-profile) to be free parameters. If the S/N is not sufficiently 
high to allow a fit, we adaptively bin up to 3 x 3” (in steps of 0.2") 
and re-fit. 

In Fig. [3] we show the velocity field of the underlying Gaus¬ 
sian emission, as well as that of the two Hi absorbers. The underly¬ 
ing emission has a strong velocity gradient (~ 700km s“^) across 
the 150 kpc in projection. At the position of the galaxy, the Lya 
is blue-shifted by ~ 150 km s“^ from the systemic, and this veloc¬ 
ity offset decreases up to ~ 1200km s“^ towards to North-West. 
The Hi absorber which has a velocity closest to the underlying 
emission, Hl-R, does not have a strong velocity gradient across 
the major axis. However, the most blue shifted absorber, Hl-B is 
blue-shifted from the underlying emission by 700 km s“^ with a 
velocity field that apparently mirrors that of the underlying Gaus¬ 
sian profile. To highlight these differences, we extract the dynamics 
through a pseudo-slit from the major kinematic axis of the halo and 
show this in Fig.|5] 

Turning to the column density maps. Fig. |4] shows that both 
of the Hi columns are remarkably uniform (although at two dis¬ 
tinct velocities), with both remaining as single systems of the full 
extent of the Lya emission. The spatial and kinematic coherence 
of the absorbers suggest they may be physically distinct from the 
background source. Using the column densities and sizes, the to¬ 
tal neutral mass of Hl-B and Hl-R are Mhi-b 10® Mq and 
Mhi-r 10® Mq. The total hydrogen mass of the absorber de¬ 


pends on the ionisation state. Photoionisation models of absorption 
systems with column densities in the range n(Hl) = 10^®“®® cm“® 
suggest average IGM neutral fractions at z = 4 of Xhi 1.5-10 
IMcOuinn et alj[201 ih . Of course, these values may represent up¬ 
per limits in an over-dense region which self-shields from UV radi¬ 
ation. Nevertheless, this suggests a total mass for the most massive 
absorber, Hl-B of ~10^° M©. 


3.3.1 Inflow or Outflow ? 


Previous observations of high-redshift radio galaxies have iden¬ 
tified very stron g (high column density) absorbers seen against 
the Lyg profiles l Hipn elein & Meisenheimerll 19931 : Ivan Oiik et al.l 


1 1 997l : iHumphrev et al.ll2007 1. The interpretation of this absorption 

is the fortuitous alignment with the intervening Lya forest absorp¬ 
tion, which may be enhanced in the biased regions around mas¬ 
sive, high-redshift galaxies. Whilst the properties of absorbing ha¬ 
los around radio galaxies have been studied in detail, most stud¬ 
ies have been limited to (either) low resolution spectra and / or 
integrated spectra, but where it has been possible to spatially re¬ 
solve the absorption, the implied covering factor is close to unity 
- on scales from 10-50 kpc - and the inferred mass of the neutral 
phase is 10^-10® M©. The absorbers also appear to be ubiqui¬ 
tous when seen against radio quiet galaxies, such as the giant “Lya 
blobs” (LABs) at z = 'i.\ which have also been studied in detail 
jMatsuda et alJ[2oTIh (we note that the other absorber in this field 
also shows a strong, spatially exten ded absorber in t he MU SE data; 
Urrutia et al. in prep). For example. iMatsuda et al.l ( l2006h find that 
the fraction of LABs with absorbers seen against the Lya emission 
is 83%. However, the nature of the absorbing gas remains poorly 
understood even when two-dimensional spectra are available. In 
particular, two-dimensional spectra are required to test whether the 
neutral gas is likely to be in-falling or out-flowing. With sufficient 
signal to noise and spectral (and spatial) resolution the spatially 
resolved dynamics of the neutral gas provided by MUSE, and its 
relation to the (underlying) emission should allow us to distinguish 
between infall and outflow, and we consider both possibilities here. 

In a scenario in which the absorber(s) around TNI 1338 arise 
due to foreground material, the giant Lya halo only provides the 
background (illuminating) screen from which the absorbers are de¬ 
coupled. In this case, the coherence of absorbers may suggest we 
are seeing a sheet (or filament) with a projected separations from 
the radio galaxy of 4 Mpc and < IMpc (assuming a Hubble con¬ 
stant at 2 = 4.1 of H( 2 ) = 423 kms“^ Mpc). The low column density 
absorber, Hl-R is presumably within the virial radius of the halo 
associated with the radio galaxy since the velocity offset places it 
within ~ 1 Mpc of the galaxy. 

Whilst this interpretation is appealing since high density 
IGM is expected to be seen around massive galaxies, there are 
a number of problems with this inflow model. Eirst, the Hi ap¬ 
pear to be enriched, although this does not preclude them be¬ 
ing unassociated with the radio galaxy, especially in an over- 
dense (biased) region. Second, the velocity structure of the un¬ 
derlying Lya emission (which has an average velocity gradient 
across the halo of dv/dR~ 4km s“^ kpc“^) is mirrored by Hl- 
B (dv/dR~ 3 km s“^ kpc“^) strongly suggesting that the material 
associated with Hl-B originated in the Lya halo (which would 
then naturally explain the similarity in the Lya and CiV line pro¬ 
files). Third, across the halo, the low density absorber, Hl-R, has 
a median Doppler 6-parameter of 6 = 170 ± 3 km s“^. We can test 
if this is “typical” of Lya abs orbers by using qua sar sight lines. 
iFechner & Reimer^ ( |2007|) and lSchave et al.l i2000l) derive column 
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density and Doppler h parameters for neutral hydrogen in the Lya 
forest at z = 2-A.5. In total, they detect ~ 1500 absorbers with 
columns between log n(Hl/cm“^) = 12-17 and & = 5-200 km s“^. 
Cutting their sample to the same median column density as seen 
in Hl-R, log n(Hl /cm“^) = 14.5, the median Doppler parameter is 
b = 30± 1 kms“^, a factor 5.5 ±0.4 lower (for the same median 
column density) as that seen in Hl-R. Thus, it appears that the low 
column density absorber seen against the halo of TNJ 1338 does 
not have the properties expected for simply being “average” IGM, 
but instead has a much higher velocity width. 

A more natural explanation for the absorbers is therefore 
that Hl-B represents the out-flowing material swept up in a wind 
shell ejected from the host galaxy and Hl-R is the low den¬ 
sity material between this shell and the host. Galactic scale out¬ 
flows from high-redshift galaxies, driven by the c ollective ef- 


fects of star-formation and supemovae winds (e.g. Pettini et al. 

2002 

; Bower et al 

2004 Erbetal. 20061: 

Swinbank et al. 

2007, 

200 s 

; [Steidel et al 

[ 20 id: [Martin et alj[2012 

: Newman et al 

2012) 

and/or AGN (e.g 

Nesvadba et al. 2007; 

Alexander et al. 

2011; 

Martin et alj|2012l: [Harrison et al.l 20121 20141; Genzel et al. 

2014) 


have been studied for some time and outflow velocities reaching 
~1000kms“^ in the ionised gas spread over tens of kilo-parsecs. 


Can we explain the properties of the absorbers in such a 
model? Numerical models can be used to trace the evolution of the 
ISM following the energy injection into an ambient medium, thus 
making predictions for th e effect of a sta rburst or AGN-driven wind 
on the ISM. For example. iKraus 3( l2005h use a hydro-dynamic code 
to model a galactic wind (due to star formation and supernova) 
on the ISM of a host galaxy and superimpose an AGN-driven jet. 
Assuming the wind starts long before the jet activity, the starburst 
first drives the wind in to a radiative bow shock, which then cools 
and is seen in Lya absorption. When the jet reaches the shell, it 
may destroys it (at least partially) so that the larger sources are 
no longer absorbed (the shell eventually cools and fragments to 
form globular clusters). In tbis model, the power required to drive 
the wind is approximately L > 5 x 10"^® erg/s ^TjlQ^K 
(Wsheii/200kms“^)^ (rsheii/25kpc). For the outflow in 
TNJ 1338, with Usheii ~1000kms“^ and rsheii ~200kpc, 
L ~ 1 X 10^® erg s“^. Although this calculation should be 
considered approximate, the luminosity required to drive the 
shell is consistent with the AGN lumino sity of TNJ 1338, 
Lagn ~2 X 10'*®ergs“^ dProuart et al.ll201^ . 

Another approach is to consider the energy of the AGN (or 
star formation) and the likely time-scales involved. Assuming the 
outflow in TNJ 1338 is a coherent structure, the mass swept up in 
the shell must be of order Moutflow ~ 10^° Mq (for a neutral frac¬ 
tion of 10%) and it will have a kinetic energy of ~10®® erg (assum¬ 
ing a spherical shell with mass ~ 10^° Mq moving at 1200 km s“^ 
with respect to the galaxy). If the shell is located a minimum of 
~ 200 kpc from the radio galaxy (this minimum being set by the 
spatial extent of the Lya halo and coherence of the absorber), 
then the minimum lifetime of the outflow is 100 Myr. Consider¬ 
ing the AGN luminosity of this source, Lagn ~ 2 x 10'*® erg s“* 
and a minimum timescale of 100 Myr, the total energy available 
is ~5 X 10®^ erg. Alternatively, considering star formation as the 
energy injection mechanism and assuming supernovae provide 
10*® erg per solar mass of stars formed, the energy provided by 
SNe will be ~ 10®® erg. Thus energetically, this star-formation or 
AGN could drive the wind if they are able to couple < 10% of their 
energy to the wind (and the energy in the wind should be consid¬ 
ered a lower limit since the galaxy-shell distance is a lower limit). 

We can also compare the mass in the absorber to the total 


ISM mass of the host galaxy. Fitting a modified black-body to the 
Herschel 250-500/rm, SCUBA 850/tm and MAMBO 1.2 mm pho - 
tometry of this galaxy dPe Breuck et al .l2004l:lDrouart et al.l2014l) . 
adopting a dust emissivity of of /3 = 1.5-2.0 ( IMagnelli et al.l 2012h 
and dust mass absorption co efficient of KsTO/im = 0.15 m® kg“^ 
dWeingartner & Drainel 1200 ih . we estimate a cold dust mass 
of Mduat '^3 X 10® Mq. For a dust-to-gas ratio of 5 dgr = 90 
dBothwell et al.ll2oT3l : ISwinbank et al.ll2014h . this suggests a H 2 
mass of Mh 2 ~3 x 10^ ® Mq, which L comparable to that esti¬ 
mated for other HzRGs dSmonts et alj|2014| ; iMilev & De Breuck 
l2008h . We caution that there is a factor 3^ x uncertainty on this 
value due to the (moderately low) signal-to-noise of the 850 fim and 
1.2 mm detections (and upper limits on the 250, 350 and 500 pim 
photometry), the choice of dust model parameters, and the uncer¬ 
tainty in the adopted dust-to-gas ratio. Nevertheless, the H 2 mass of 
the ISM in the radio galaxy is similar to that estimated for the total 
mass of the two structures we see in absorption against the halo. 

Will the gas escape the halo potential? If we consider a 
-^10** Mq galaxy inside a dark halo with mass ~10*®-10** with 
NFW dNavarro et alJll9^ density profile, the gas would require 
a velocity greate r than ^ 1000kms~* to escape the halo (see also 
the calculation in iHarrison et al.l2012h . This is similar to the veloc¬ 
ity seen in Hl-B. However, some models have suggested that even 
massive outflows may stall in the galaxy halo (or in this case, the 
“group” halo), and re-collapse and cool at later times (along with 
new fuel supplies), resulting i n re-ignition of sta r formation and fur- 
ther black-hol e grow th (e.g. iLagos et'^l2008l : iGabor et'^l201 ll : 
iHopkins et al.l 2013h. especially in g roup/cluster (or proto-cluster) 
environments (IMcCarthv et alj201 ih . 

In summary, while it is not possible to constrain definitively 
whether the AGN or star-formation have driven the outflow, or 
the ultimate fate of out-flowing gas, it appears that the Hi ab¬ 
sorbers are associated with activity occurring within the galaxy, 
and with outflow properties (high velocity, spatially extended, and 
mass loaded) in broad agreement with the requirements of mod¬ 
els that invoke AGN-driven outflows to r egulate star formation 
and black-hole y owth in luminous AGN dpi Matteo et alj[2005l : 
iHopkins & Elvia[20ld : IPebuhr et al.ll2012h . Even if the outflow 
does not escape the galaxy halos, it may have sufficient energy 
to heat the halo and co ntrol the level of coo l ing in massive ha- 
los at later times (e. g. oiurazov et alJ 1200^ : Iciabor et al.l I 201 I] : 
[McCarthy et al.ll201 ll-lBower et al.ll2012l) . 


4 CONCLUSIONS 

We have used the MUSE integral field spectrograph to study the 
giant Lya has around the high-redshift radio galaxy TNJ 1338 at 
2 = 4.11. Our main observations reveal the following: 

• Our observations map the intensity distribution and dynamics 
of the Lya halo which is extended over a region approximately 
150 X 80 kpc. In comparison, the radio jets are confined to the cen¬ 
tral 10 kpc of the halo, although the alignment of the extended 
Lya halo and radio jets may be indicative of previous AGN activ¬ 
ity (which ionised the extended halo and which is now seen in Lya 
emission as the gas cools). 

• The spatial extent of the Lya halo in TNJJ1338 is one of the 
largest seen in high-redshift radio galaxies. We clearly detect the 
Civ emission in the extended halo suggesting that the gas in the 
extended is enriched. 

• Across the entire halo, the Lya emission line appears asymmet¬ 
ric. The line profile is best-fit by two H-I absorbers superimposed 
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Figure 5. The velocity profile of the underlying Gaussian emission and the 
two absorbers extracted across the major axis of the Lyo: halo (as shown in 
Fig.ig. The solid points and thick eiTor bars denote the central velocity and 
the en'or on the central velocity, whilst the thinner errors at the top of the 
panel show the average width of the Gaussian (cr) or Doppler 6-parameter 
(absorbers) in each bin. The origin of the velocity and spatial axes is fixed 
on the radio galaxy, demonstrating that these motions do not reflect bulk 
rotation o f gas around the radio galaxy if it lies at the bottom of the potential 
well (c.f. Ivan Oiik et al .I 1997 I) . 


on an underlying Gaussian emission profile. One of the striking 
results from our MUSE observations is the coherence and spatial 
extent of both absorbers. Indeed, both absorbers have a covering 
factor of unity. The first absorber, Hl-B appears blue-shifted from 
the (underlying) Gaussian emission line profile by ~ —700km s“^ 
and has a column density of n(Hl)~ 10^®'^ cm“^. Since the veloc¬ 
ity structure of this absorber follows that of the underlying emis¬ 
sion, this material is most likely associated with an out-flowing 
shell of material, and we estimate the mass in the shell must be 
of order ~10'^° M©. By comparison to the AGN luminosity or star- 
formation rate of the galaxy, we find that both could trivially drive 
the wind if they are able to couple < 10% of their energy to the 
wind. The second absorber, Hl-R is much closer to the underlying 
halo emission. Aw ~ —200kms“^ but has a much lower column 
density, n(Hl)~ cm“^. The Doppler b parameter for this ab¬ 
sorber is ~6x higher than expected for the low-density IGM (at 
the same column density), suggesting that this is not just a typical 
intergalactic absorber, but instead that this material represents the 
low-density IGM between the out-flowing shell and the Lya halo. 


We show that the energetics of the star-formation or AGN could 
drive the outflow if it is able to couple a small fraction of tbe en¬ 
ergy to the gas. 

• The underlying Gaussian emission line profile shows a velocity 
gradient of '~700kms“^ across the 150 kpc in projection. At the 
position of the galaxy, the Lya is blue-shifted by ~ 150kms~^ 
from the systemic. The velocity dispersion profile of the underlying 
emission (which has an average value of (t = 570± 13kms“^) is 
peaked at the position of the radio galaxy and the radio-jets, which 
may reflect the interaction of the radio plasma and the gas. 


• We identify two other 2 = 4.1 emitters in the field of the ra¬ 
dio galaxy which are presumably associated with the group/cluster 
within this over-density. We also serendipitously identify a new gi¬ 
ant Lya halo at 2 = 3.31 approximately ~ 30” to the south west of 
TNJ J1338 which appears to be associated with a number of other 
star-forming galaxies and AGN at this redshift. 


Finally, we can consider how the galaxy and outflow may 
evolve towards the present day. The out-flowing shell will be both 
thermally and gravitationally unstable, and will eventually frag¬ 
ment. At this point, the covering factor will decrease and the ab¬ 
sorption will vanish (which would explain why no absorbers are 
seen in the halos of radio gala xies with jet activity > 50 kpc scales; 
[Rottgering & Pentericci| l999|). The Jeans mass in the shell will be 
of order 10® Mq ( Krausell2002li . Hence the shell will fragment and 
form stars in globular clusters of this mass. Of course, at this point 
the shell should become visible in the rest-frame UV / optical as 
young stars form (which will also increase the ionisation of the 
remaining parts of the shell). Assuming an efficiency of 10%, 
this shell will contain a ~ 1000 globular clusters of this mass - in 
reasonable agreement wit h the number of cl usters found in nearby 
brightest cluster galaxies dHarris et al .Il998h which show excess of 
globular clusters systems compared to non-BCGs of comparable 
luminosity. 

In summary, the most successful models of galaxy forma¬ 
tion increasingly require outflows of gas to expel gas from galax¬ 
ies at early times. In the most massive galaxies, most of the en¬ 
ergy injection is believed to arise from the central AGN, making 
high-redshift radio galaxies, which have many of the character¬ 
istics expected for the progenitors of todays massive Ellipticals, 
strong candidates to empirically constrain the properties of out¬ 
flows which can be used to test the models. The discovery of high- 
velocity (~ 1000kms“^), halo-scale (i.e. 150kpc) and mass 

loaded winds in the vicinity of the central radio source are broadly 
in agreement with the requirements of models that invoke AGN- 
driven outflows to regulate star formation and black-hole growth 
in luminous AGN dOi Matteo et al.ll200^ ; iHonkins & Elvigl201(]| ; 
IPebuhr et al.ll20l3) . 
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